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We performed a detailed study of the 2007 outburst of the 352 s pulsar SAX J2103. 5+4545, a 
Be/X-ray transient observed by INTEGRAL, to study its spectral and temporal properties during 
the evolution of the outburst. SAX J2 103. 5+4545 was observed with IBIS/ISGRI from 25 to 
27 April 2007 and from 6 to 8 May 2007. The 20 - 100 keV spectrum is well described by 
a bremsstrahlung model with a temperature kT ^ 24 keV. The pulse profiles are variable with 
time and energy. A pulse period derivative of Pspin-up — 

(-3.4 ±0.4) X lO"'^ s s-^ has been 
observed during the outburst. Instead, a spin-down of Pspin-down = (5-5 ±0.4) x 10^^ s 5^' is 
observed between the 2007 outburst reported here and the previous one occurred in December 
2004. This is the largest spin-down measured for SAX J2 103. 5+4545 since its discovery. We 
estimate a neutron star magnetic field in the range (1.6 — 3) x lO'^ G using the Ghosh & Lamb 
torque model. 
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5.4210x10* 5.4215x10' 5.4220x10' 5.4225x10' 5.4230x10' 
MJD 

Figure 1: RXTE/ASM light-curve of SAX J2103.5+4545 in the 2-10 keV energy range (54210 - 54230 
MJD, April 20 - May 10 2007). The arrows indicate the time interval covered by the INTEGRAL observa- 
tions. 



1. Introduction 

SAX J2103. 5-1-4545 is a transient Be/X-ray binary with an orbital period of 12.67 days [Q] 
hosting a pulsar with a pulse period of 352.2 s [IT]]. It was discovered during an outburst in February 



1997 [11], and since then numerous type I and type II outbursts have been observed \U\. The optical 



counterpart is a BOVe star with highly variable Ha emission, located at a distance of 6.5 ± 0.9 kpc 



[]14[]. The position of SAX J2103. 5-1-4545 in the spin period versus orbital period diagram deviates 
from the correlation of HMXB with Be star and falls in the region typical of wind-fed HMXB with 
OB supergiant companions. The X-ray energy spectrum in the 5-200 keV band is described by 
a power-law with exponential cut-off, with photon index F = 1.53, cut-off energy Ec = 19 keV 



and e-folding energy Ef = 32 keV. No cyclotron lines have been observed [|16|]. A soft blackbody 
component with temperature kT ~ 1.9 keV was observed with XMM-Newton [12], and an iron 
emission line at ~ 6.42 keV was observed withRXTE [||]. SAX J2103. 5-1-4545 has shown numerous 
spin-up phases during its outbursts. A correlation between X-ray flux and spin-up rate has been 
found [0] which can be explained with the formation of an accretion disk. 



2. Observations and data analysis 

SAX J2103. 5-1-4545 underwent an outburst in April-May 2007 [Q]. This outburst was observed 
by the IBIS/ISGRI instrument on 25-27 April 2007 (54215.39 - 54217.02 MJD) and on 6-8 May 
2007 (54226.71-54228.98 MJD). 

Here we report the data analysis of the source field observed with IBIS-ISGRI. The Off-line 
Scientific Analysis (OSA 6.0) software was used. We have divided this observations in 2 groups 
(Table [I|) which are marked in Figure where the 2-10 keV light-curve of the 2007 outburst 
observed by RXTE/ASM is reported. 

We extracted two spectra of SAX J2 103.5-1-4545 in the energy range 22 — 100 keV, in the time 
intervals 25-27 April 2007 and 6-8 May 2007 (Table [T]). A systematic uncertainty of 1% was added 
quadratically to the statistical errors to account for uncertainties in the instrumental response. 
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Obs. 


Time interval 


Start Time 


End Time 


Exposure time 






(MJD) 


(MJD) 


(ks) 


1 


25-27 April 2007 


54215.39 


54217.02 


129 


2 


6-8 May 2007 


54226.71 


54228.98 


169 



Table 1: The two groups of IBIS/ISGRI observations of SAX J2103.5+4545 used in this work. 



Time interval 


kT 


Xv 


Flux (20 - 40 keV) 


Flux (20 - 100 keV) 




(keV) 


for 13 (d.o.f.) 


10^^ erg cm^2 s^' 


10^^ erg cm^2 s^^ 


25-27 April 2007 


23.4+^;^ 


1.06 


, 1 1 +0.04 
^■^^ -0.03 


1.68 _Q 


6-8 May 2007 


24+^ 


1.29 


0.32 ±0.03 


0.50 ±0.06 



Table 2: Best fit parameters of the two average spectra of SAX J2 103 .5+4545 with a bremsstrahlung model. 

We performed timing analysis in the range 15 — 60 keV, where the best signal to noise ratio is 
obtained. We corrected the arrival times to the Solar System barycenter and for Doppler delay due 
to the pulsar orbital motion. The correction of the arrival times can be written as: 

3 1 1 

5t =xsml — -.xe sin CO ± -;ice cos CO sin (2/ ) — sin CO cos(2Z) (2.1) 

where x = sin(/)/c is the light-travel time for the projected semi-major axis Qx, i is the inclination 
angle, e is the eccentricity and CO is the longitude of periastron [^. Z is the mean orbital longitude 
at time t which is given by: 

where Tj^^2 is the epoch when the mean orbital longitude is equal to 7r/2 and Port is the orbital 
period. To correct the arrival times, we have assumed the most recent orbital parameters obtained 
by Camero Arranz et al. (2007) [^. An epoch folding technique was used to measure the pulse 
period |[l|]. 

3. Results 

3.1 Spectral analysis 

We have fitted the 22 — 100 keV spectra with several simple models and found that the best 
X'^ value is given by a bremsstrahlung model with a temperature kT 24 keV for both spectra (see 
Figure ^ for energy spectra, and Table ^ for parameters). 

The average flux (20-100 keV) obtained from the spectrum in the time interval 25-27 April 
2007 is = (1 .68 m ) x 10"*^ erg cm^^ , while on 6-8 May 2007 the flux is = (5.0+° og ) x 
10^'" erg cm^2 s^\ translating into a luminosity = (8.5 ±2.5) x 10^^ erg s^' (25-27 April 
2007) and = (2.5 ± 0.7) x 10^^^ erg s"' (6-8 May 2007) for a distance d = 6.5± 0.9 kpc. These 
luminosities are similar to those measured during previous outbursts (|jT6|l, [|I]]). 
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Figure 2: ISGRI SAX J2103.5+4545 spectra in the energy range 22-100 keV, on 25-27 April 2007 (Top) 
and on 6-8 May 2007 (Bottom). Left panels: counts spectra, best-fit bremsstrahlung spectra and residuals 
(in units of standard deviations); Right panels: photon spectra together with the best-fit model. 



Name Time interval Pspm (s) 

A 25-27 April 2007 352.725 ±0.004 
P2 6-8 May 2007 352.382 ±0.004 



Table 3: Pulse periods of SAX J2103.5h-4545. 



3.2 Timing analysis 

We measured the pulse period and the associated uncertainty for each of the two data sets 
reported in Table |l| by means of an epoch folding technique, obtaining the values reported in Table 
^ The average spin-up rate between Pi and P2 is Pspin-up = (—3.4 ±0.4) x 10~^ s s^' The average 
flux during this spin-up is = 3.6 x 10^^ erg cm^^ s^^ in the energy range 1 — 200 keV. 

SAX J2 103.5-1-4545 has shown a pulse period decrease since 1997, occasionally interrupted 
by spin-down intervals, corresponding to low accretion rates, ([p^, [jl]]). In particular, the pulse 
period history (Figure ^ shows a spin-down between December 2004 (53340 ± 10 MJD) where 
PBaykai 2007 = 352.31 ±0.01 s and the last outburst (April - May 2007, see Table [I]), when 
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Figure 3: Pulse period history of SAX J2103. 5+4545. The pulse period values are taken by BeppoSAX Ql 1| 
INTEGRAL and RXTE 



we measured Pi = 352.725 ± 0.004 s. Between these two outbursts SAX J2103.5+4545 was in 
a low luminosity state. We measured the upper limit flux of SAX J2103. 5+4545 during this low 
luminosity state from December 2004 to April 2007 using RXTE/ASM data. We found F <2x 
10^'^ erg cm^^ s^^ in the energy range 2 — 10 keV. The average pulse period derivative from 
December 2004 to April 2007 is Pspin-down = (5.5 ±0.4) x 10"'^ s s-K 

For each of the 2 data sets reported in Table |l], we have obtained the pulse profiles for the 
energy bands 15-60 keV, 15^0 keV and 40-60 keV. The pulse profiles are binned into 20 phase 
bins (zero phase at 51000 MJD; see Figure 0). We have found a temporal variability for the pulse 
profiles, as already observed during previous outbursts of this source ([^, [|l6|]). 



4. Discussion and Conclusions 

SAX J2103. 5+4545 is a Be/X-ray binary with transient X-ray emission, fed by the wind of the 
Be companion. In previous outbursts a correlation between X-ray flux and spin-up rate was found 
which can be explained with the formation of an accretion disk [^. 

When a neutron star is fed by an accretion disk, the material which flows towards the neu- 
tron star magnetic poles produces X-ray luminosity and at the same time a torque on the rotating 
neutron stars mediated by its magnetosphere. The angular momentum of the infalling matter is 



transferred to the pulsar at the magnetospheric boundary [IC]. The Ghosh and Lamb model as- 
sumes a magnetically-threaded disk allowing not only spin-up, but also spin-down, which occurs 
when the magnetic and corotation radii are comparable and the luminosity is low. Ghosh and Lamb 
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Figure 4: Evolution of SAX J2103. 5+4545 pulse profiles for different time intervals and for two energy 
bands. From top to bottom: 15 - 40 keV, 40 - 60 keV. From left to right: Obs. 1, Obs. 2 (see Table 3). 



predicted the following correlation between the pulse period derivative and the X-ray luminosity: 



P = -5.0 X lO-VdM<o.)Rri.s' iPLff 



-3/7 



s yr 



(4.1) 



where is the neutron star magnetic dipole moment in the disk plane {pL = BE?) in units of 
lO^o G cnv', /?6 is the stellar radius in units of 10^ cm, Mns is the mass of the neutron star in 
grams, 745 is the moment of inertia in units of 10^^ g crn^, P is the pulse period, L37 is the accretion 
luminosity in units of 10^'' erg s^^. We have considered the dimensionless torque obtained by 



Ghosh & Lamb (1979) [IC] 



n{co,) ^ 1.39 



1 - [4.03(1 - w,)'' "3 - 0.878] 



1 -O), 

where (Os is called fastness parameter and is defined as: 



3/2 



(4.2) 



(4.3) 



where Q.s is the angular velocity of the star, Q.K{ro) = {GMf^s/^IY^'^ is the Keplerian angular 
velocity at ro, the radius which divides the outer part of the accretion disk, where the velocity is 
keplerian, from the inner part, where the plasma corotates with the star [10|; Rc = {GMns/^^Y^^ 



is the corotation radius of the star. Therefore, during an outburst, a variation in the accretion 
rate can modify the pulse period derivative and at the same time the X-ray luminosity. In this 
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Figure 5: Comparison between the observed spin-up and spin-down rate and the theoretical relation between 
P and the source luminosity. The box labeled with Pspin-up indicates the spin-up rate measured during the 
time interval 25 April - 8 May 2007, the box labeled with Pspin-down indicates the spin-down rate measured 
during the time interval December 2004 - April 2007. The fit with theoretical curves is obtained for a 
distance ~ 4.8 kpc and a magnetic field B= 1.6x lO'^G (soUd line) and B = 3xlO'^G (dashed hne). The 
uncertainties of Pspm-up, Pspin-dmvn and are given by the pulse period derivative uncertainties and by the 
uncertainty of the source distance. 



framework it is possible to determine the magnetic field of SAX J2103.5-I-4545. Using equations 
( P~T| ) to (13), we assumed Myvs = 1.4 Mq, Rns = 10 km> an average pulse period P = 352.5 s. 
Fx = 3.6 X 10^^ erg cm^^ s^^ as average flux during the spin-up, and the upper limit flux during 
the spin-down. We found a magnetic field of (1.6 — 3) x 10^^ G for a distance of --^ 4.8 kpc. This 
distance is consistent with what found by Reig et al. (2004), if the error in the absolute magnitude 
found by Vacca et al. (1996) JT/I ] is included. These results are also consistent with what previously 
obtained by Baykal et al. (2007) [0]. In Fi gure (||) we compare the spin-up and spin-down rate 
observed with INTEGRAL with the theoretical relation between P and the source luminosity for 
3 different magnetic fields of the neutron star. Vertical lines mark the source luminosities for 
d = 4.8 kpc and d = 9.2 kpc. 

During the low luminosity state (December 2004 - April 2007) the upper-limit to the luminos- 
ity is Ljf < 5 X 10^^ erg s^^ (2 — 10 keV, d = 4.5 kpc). This low luminosity cannot be explained with 
the high-density (M = 10^^ Mq yr^^) wind of the circumstellar disk around the Be star. From op- 
tical observation during low luminosity states, SAX J2103. 5-1-4545 shows Ha in absorption, which 
is interpreted as the loss of the circumstellar disk (Q, [||]). In this case, SAX J2103. 5-1-4545 
could accrete from the low-density wind of the polar regions. The mass loss rate out of circum- 
stellar disk is 10^^ — 1O^^M0 yr^^, the velocity law is v = Voo(l —RBe/r)^, with terminal velocity 
v<» ^ 600-1800 kms-\^^ \ ([|^], [0]). Assuming the SAX J2103.5-i4545 parameters found 
by Camero Arranz et al. 2007 [H], and assuming Mbc = 20 Mq, R^s = 10 km, Mf^s = 1-4 Mq and 
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Rbs = 8 Rq ([14], [15]), we find that the expected luminosity in the Bondi-Hoyle accretion theory 
([§, [0) is L, < 10^4 

erg s ^ If we consider the pulsar magnetic field of 3 x 10^'' G, calculated 
above, a wind mass loss rate of 8 x 10^^ Mq yr^^, and Vco ~ 600 km s^^, the centrifugal inibition 
of accretion produces an expected luminosity = ^^^M < lO^^erg s^^, where R,„ is the mag- 
netospheric radius [0]. These expected luminosities agree with the measured upper-limit, and with 
the observed spin-down rate Pspin-down = (5.5 ±0.4) x 10^^ s s^^. 
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